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C$ ■ Abstract 



Using Chandra X-ray, Spitzer mid-IR, and 1.5 GHz radio data, we examine the 
spatial structure of SNR 3C 391. The X-ray surface brightness is generally anti- 
correlative with the IR and radio brightness. The multiband data clearly exhibit 
a heart-shaped morphology and show the multi-shell structure of the remnant. A 
thin brace-like shell on the south detected at 24 /xm is projected outside the radio 
border and confines the southern faint X-ray emission. The leading 24 /mi knot on 
the SE boundary appears to be partly surrounded by soft X-ray emitting gas. The 
mid-IR emission is dominated by the contribution of the shocked dust grains, which 
t^- ', may have been partly destroyed by sputtering. 
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1 Introduction 



Galactic molecular clouds, containing about half the mass of the interstellar 
medium (ISM), have a complex structure, and stars form out of the densest 
regions therein. Massive stars have often not moved far from their birthsites 
by the time they explode as a core-collapse supernova, and therefore the in- 
teraction of supernova remnants (SNRs) with molecular clouds can plausibly 
occur. The number of SNRs with convincing evidence for interaction with 
ambient molecular clouds is about 20, and a large fraction of them are mixed- 
morphology or thermal composite SNRs, with interior-filled thermal X-ray 
emission surrounded by a radio shell (Yusef-Zadeh et al., 2003). 

When SNRs interact with inhomogeneous molecular gas with a large density 
variation, the morphology and structure is not expected to be regular or sym- 
metric. Features such as gas clumps, incomplete arcs, blowout shells, thin or 
thick filaments, etc., can typically be observed in radio, infrared (IR), optical, 
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or X-rays. SNR 3C 391 is one of the prototype mixed-morphology or thermal 
composite SNRs, confirmed to be interacting with molecular clouds to the 
northwest (Rho & Peter, 1998; Wilner et al, 1998; Reach & Rho, 1999; Chen & 
Slane, 2001). Due to its remarkable radio shell elongated from northwest (NW) 
to southeast (SE), it is suggested to have broken out of a dense region into an 
adjacent region of low density (Reynolds & Moffett, 1993). The bright clumpy 
thermal X-ray emission arises from the remnant interior, and the rim is X-ray 
faint (Chen & Slane, 2001; Chen et al., 2004, hereafter CSSW04). The Spitzer 
Space Telescope (SST) (Werner et al., 2004) Infrared Array Camera (IRAC) 
near-IR observations show a bright incomplete shell correlated well with the 
NW radio structure which is tangent to a giant molecular cloud (Reach et 
al., 2005; Lee, 2005). These structure patterns at different wavelengths show 
complex physical conditions in SNR 3C 391, and reflect inhomogeneities in 
the ISM. 

Here we will show as overall view of SNR 3C 391 combining the radio (1.5 
GHz), mid-IR (SST 24 & 70 /im), and X-ray (Chandra) observations, and 
examine the spatial structures in different wavelengths. 



2 Multi-band observation of 3C 391 

2. 1 The Radio and X-ray Morphology 

In the radio images (Reynolds & Moffett, 1993; Moffett & Reynolds, 1994), 
3C 391 has a bright partial shell along the north and NW boundary. The 
NW part of the remnant appears to be a broken bubble, with the brightest 
radio bar along the western border and a faint blowout extending to the SE. 
Two 1720 MHz OH maser spots are detected along the radio shell (Frail et 
al., 1996). Reynolds & Moffett (1993) suggest that the remnant may have 
resulted from a supernova explosion in the interior of the molecular cloud and 
is expanding into a region of greatly varying density. 

The Einstein (Wang k Seward, 1984) and ROSAT (Rho & Peter, 1996) X- 
ray observations show a strong soft X-ray emission arising from the southeast 
region of the remnant. Both the ASCA and Chandra observations show a 
SE-NW elongated broad-band X-ray morphology, similar to the radio image 
(Chen & Slane, 2001; CSSW04). The X-ray emission peaks southeast of the 
geometric center of the radio pattern. The X-ray emission from the interior 
to the NW bubble suffers higher absorption than that from the SE blowout 
part. The Chandra observation reveals rich small-scale structures, such as 
knots and arcs, as well as faint, diffuse gas that appears to expand out of the 
southwestern radio boundary (CSSW04). In the southwest (SW), the narrow 



2 



band analysis of the S and Si X-ray lines unveils a very faint jet-like protrusion 
projected outside the radio border (Su & Chen, 2005). 

To better show the entire spatial X-ray emitting remnant, we rebin the cleaned 
level 2 Chandra data (ObsID 2786; CSSW04) using an adaptive mesh to in- 
clude at least 9 counts (0.3-7 keV) per bin. The contours of the rebinned 
image are plotted in Fig. la. In this figure, a broad band X-ray protrusion 
(corresponding to the narrow band one) is seen in the SW. 



2.2 Mid-IR images 

The mid-IR 24 /im and 70 /im observations used here were carried out as part 
of the 24 and 70 Micron Survey of the Inner Galactic Disk Program (PID: 
20597, PI: Sean Carey) with the Multiband Imaging Photometer (MIPS) 
(Rieke et al., 2004). The data are obtained from Spitzer archive. The raw 
24 /im image of 3C 391 from the Post Basic Calibrated Data is shown in 
Fig. lb, which displays the remnant extent which is similar to that seen in 
radio and X-rays. We sum 24296 pixels ( ~ 3.7 arcmin radius) to get the 
total surface brightness and select faint regions outside the SNR as the back- 
ground. Subtracting the mean surrounding IR background ~ 60 ± 3 MJy/sr, 
we obtain a mean surface brightness of 3C 391 ~ 11.5 ± 3.0 MJy/sr or a 
flux ~ 37.9 ± 9.9 Jy. After a manual subtraction of nearby emission fea- 
tures, the remaining 24 /im emission is superposed on the X-ray contours in 
Fig. la. The manually subtracted 24 /im and 70 /im emission maps are over- 
laid with the 1.5 GHz radio contours in Fig. 2. Tricolor image with 1.5 GHz 
radio emission in red, 24 /an IR emission in green, and 0.3-7.0 keV X-rays in 
blue is shown in Fig. 3. In the 70 /im mid-IR image there are some streaks 
mixed with the IR emission and it is very difficult to determine the back- 
ground intensity or the pure mid-IR intensity associated with SNR. We use 
the GeRT packa ge (|http://ssc.spitzer.caltech.edu/mips/gert[) and the software 
tool MOPEX (http:/ /ssc. spitzer. caltech.edu/postbcd/download-mopex.html) 
to creat a new mosaic of 70 /im. After running MOPEX on the filtered Ba- 
sic Calibrated Data (BCD) images produced by the data processing pipeline, 
some arc-like structures of 3C 391 are visible in 70 /im (Fig. 2b). 

The two mid-IR images (Fig. 2) display several thick arc-like structures along 
the SW, NW, and northeast (NE) boundarier. A very faint 24 /im brace- 
symbol-like thin arc structure in the south is projected outside the radio border 
and goes from west to south along the edge of the faint X-ray emission. 

In the 24 /im emission map (Fig. la, lb), there are two bright knots at 
(18 h 49 m 39 s .9, -00°56'40") and (18 h 49 m 39 s .4, -00°58'36") on the east and south- 
east border, respectively. They seem to correspond to X-ray knots (Fig. 3). 
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Fig. 1. (a) Superposition of Chandra X-ray (0.3-7.0 keV) contours of 3C 391 on 
the MIPS 24 /im image. The X-ray contours are between the 5% maximum and 
the 90% maximum with a square-root scale. The IR intensity range is from 62 MJy 
sr _1 to 80 MJy sr -1 , with a linear transfer function and the point-like IR sources 
away from the SNR have been removed, (b) The SST MIPS 24 /xm raw image of 
SNR 3C 391. The IR intensity range is from 58 MJy sr" 1 to 92 MJy sr" 1 , with 
a linear transfer function but the background and the point-like sources have not 
been removed. 



The east knot is the 24 fim counterpart of IR source IRAS 18470—0100 
(CSSW04), with improved positional resolution. It can also be seen that a 
line of point-like IR sourses are evenly spaced extending from the remnant 
center to the bright SE knot. These five sources have 2MASS near-IR coun- 
terparts with accurate positions: (18 h 49 m 27M7, -00°56'44".l), (18 h 49 m 29 s .34, 
-00°57'29".9), (18 h 49 m 31 s .43, -00°58'02".5), (18 h 49 m 35M5, -00°58'26".l), 
and (18 h 49 m 39 s .50, -00°58'36".0). The IR source at the geometric center has 
an X-ray counterpart (CXJ 184927.0-005640; CSSW04) given ~ 4" spatial res- 
olution of 2MASS. The nature of these sources is interesting but unknown yet. 
If the central one suffers an interstellar absorption similar to that of 3C 391 
(Nh ~ 3 x 10 22 cm -2 ), assuming its X-ray spectrum as a thermal plasma 
with a temperature of 0.6 keV (for soft case) and 2 keV (for hard case), re- 
spectively, typical for normal early-type stars, we estimate the 0.3 — 8 keV 
luminosity of the X-ray source as ~ 4.3 x 10 32 ergs s _1 (for soft case) and 
~ 3.3 x 10 33 ergs s -1 (for hard case) based on the X-ray count rates (Ta- 
ble 1 in CSSW04). It indicates that, on this assumption, the central source 
would most probably represent an individual early-type star or even a massive 
colliding stellar wind binary. 
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Fig. 2. (a) The SST MIPS 24 /im image of SNR 3C 391, overlaid with the 1.5 
GHz radio contours. The IR intensity range is from 62 MJy sr _1 to 80 MJy sr _1 , 
with a linear transfer function and the point-like IR sources away the SNR are 
removed. The overlaid 1.5 GHz radio contours are at 1.5, 3.19, 8.25, 16.69, 43.68, 
and 62.24 x 1(T 3 Jy/beam (from Moffett & Reynolds 1994). (b) The MIPS 70 fim 
image of SNR 3C 391 overlaid with the same radio contours as used in (a). The 
radio character of the image is same as in (a). The dashed lines indicate pronounced 
IR arc/shell- like structures. 

3 Discussion 



3. 1 Spatial comparison among X-ray, mid-IR, and radio emissions 

Using the available radio, IR, and X-ray data, we can investigate the struc- 
ture of SNR 3C 391. Seen in the tricolor image (Fig. 3), the X-ray surface 
brightness is generally anti-correlative with the IR and radio brightness, or 
anti-correlative with the environmental density because of intervening extinc- 
tion. In the NW half of the remnant, two distinct thick arc-like structures in 
both 24 and 70 /im coincide grossly with the radio shell (Fig. 2 and Fig. 3). 
The westernmost bar-like radio peak and the "broad molecular line" (BML) 
region are located along the west-southwest IR shell, which represents a dense 
shocked region. The two 1720 MHz OH maser spots (Frail et al., 1996) are 
harbored in the NE and SW thick radio/mid-IR shells. On the NE boundary, 
a thick mid-IR arc is coincident with the radio shell, too. These arcs/shells are 
also discernible in the SST 5.8 //m near-IR images (Reach et al., 2005; Lee, 
2005). 

Apart from the above general comparison, we note two features of particular 
interest. 
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Fig. 3. Tricolor image of SNR 3C 391. The 1.5 GHz radio is coded in red, the SST 
24 jim IR in green, and the Chandra X-ray in blue. The intensity map is square-root 
scaled in each energy band. 

1. The knot on the southeastern border (18 h 49 m 39 s .4, -00°58'36") (the SE 
knot in Fig. 2a) appears to be partly surrounded by soft X-rays (Fig. 4). Two 
possible mechanisms may explain this small-scale phenomenon. First, after 
the blast strikes a small dense cloud, a transmitted shock heats the cloud to 
emit 24/im emission, while a bow-like reflected shock stands around the cloud 
(McKee & Cowie, 1975). Second, after the cloud is shocked by the blast wave, 
materials are evaporated from the cloud and form a soft X-ray emitting gas 
layer (Cowie & McKee, 1977). 

2. The 24 /xm brace-like thin shell (Fig. 2a) appears to be projected outside 
the long southern radio border of the remnant (Fig. 2a and Fig. 3). It has not 
been seen before at any wavelengths. However, it confines very well the faint 
X-ray emitting emission that fades out from the remnant interior and seems to 
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Fig. 4. The southeastern boundary region. The 24 /jm IR emission is coded in 
green, and the Chandra soft (0.3-2.6 keV) X-ray emission in blue. 

extend beyond the SW radio border (CSSW04) (Fig. la and Fig. 3). Moreover, 
the position of the cusp of the "brace" (Fig. 2a) is coincident with a jet-like 
protrusion seen in the narrow band X-ray maps for S and Si (Su & Chen, 
2005). Such positional consistency between the 24 /xm and X-ray emission 
suggests that the thin brace-like IR shell is associated with the remnant. Thus, 
the 3C 391 SNR has a larger volume than previously recognized and takes a 
heart-like shape with combination of the X-ray, IR, and radio emission (Fig. 2a 
and Fig. 3). The multi-shells, especially the newly detected southern thin 
shell, indicate that the molecular clouds around this remnant are distributed 
more complexly than previously assumed. The NW clouds have a non-uniform 
structure and there is a low density region in the south, into which a part of 
the supernova blast wave propagates. 



3.2 The 24 [im mid-IR emission of the multi-shells 

The 24 (Urn emission more distinctly demonstrates shell structures of SNR 3C 391 
than other IR wavelength emission. Similar phenomena occur in some other 
SNRs (Stanimirovic et al, 2005; Morris et al., 2006; Borkowski et al., 2006). 
Hence the 24 fim band seems to be an effective window for SNR observation. 
The near-IR emission from various portions of the 3C 391 SNR has been sug- 
gested to be dominated by atomic fine-structure lines and synchrotron emis- 
sion (from electrons in the moderate-density molecular regions) or molecular 
lines and ionic lines (from the high density molecular gas) (Reach et al., 2005). 
Here we suggest that the mid-IR emission of 3C 391 mainly arises from dust 
grains. 
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The origin of the 24 mid-IR emission is not clear. We explored possible 
mechanisms of synchrotron, atomic lines, molecular lines, or IR emission from 
shock-heated dust grains. Let us compare the contribution of each mechanism 
to the 24 /mi emission of 3C 391. (a) Synchrotron. Using the integrated flux 
at 330 MHz of 41.0 Jy and the spectral index - -0.49±0.1 (Brogan et al. 
2005), we estimate the flux by extrapolation to be less than 0.3 Jy (1% total) 
at 24 /mi. Certainly, the mid-IR emission of the brace-like thin shell is not 
synchrotron, otherwise it should have been discernible at 70 jum. (b) Ionic 
lines. The SST MIPS band at 24 /mi contains ionic lines of [Fell] at 25.988 /mi 
and [OIV] at 25.890 /im (Reach & Rho, 2000). In the BML region and the 
radio bar-like peak, the [Fell] plus [OIV] lines amount to no more than 40% 
and 30% of the local MIPS IR 24 /im fluxes, respectively (Reach & Rho, 2000). 
We assume that similar fractions hold for the entire remnant but note that 
this may overestimate the fraction because the two regions are the brightest 
in the ionic lines, (c) Molecular lines. Some rotational molecular lines such 
as from OH and H2O may contribute to the mid-IR emission, but generally 
mid-IR from molecular lines is weak compared with the strong continuum 
emission of dust grains (Gorti & Hollenbach, 2004). Also, the ISO observations 
of this remnant do not show strong molecular lines at 24 /im (e.g., Reach 
& Rho, 2000). (d) Dust emission. Since the 24 /im emission of 3C 391 is 
not dominated by any of the above three mechanisms, we suggest that it 
predominantly (> 60%) originates from dust grains. The concentration of the 
24 /im emission on the remnant shells also favors that the dust in SNR 3C 391 
is part of the swept-up, shock-heated ISM. This physics of explanation is 
similar to that for the mid-IR emission of SNR N132D (Tappe et al, 2006). 

The blackbody temperature of the dust grains in SNR 3C 391 was estimated 
as 140 K for 12-24 /tm and 52 K for 24-100 /mi from the IRAS observation 
(Arendt 1989). The total mass M d of the hot dust associated with the SNR 
at temperature Td is given by: 

M _ m 
Md «(A)B„(A, T d ) ' [L) 



where F V (X) is the observed 24 /iin flux ( ^ 23 Jy), d ~8.0 kpc the distance to 
the remnant (Chen & Slane 2001), k x = 2.5(A/450/im)~ 2 - cm 2 g" 1 the dust 
mass absorption coefficient (Draine & Lee 1984), and B U (X, Td) the Planck 
function evaluated at the dust temperature. Thus we get a lower limit on the 
dust mass M d > 2 x 10" 4 M for T d ~ 140 K and M d > 0.3 M for T d ~ 52 K. 
With the mass of the X-ray emitting gas of the remnant M g ~ 100 M 
(CSSW04), the dust-to-gas mass ratio is then >2 x 10~ 6 -3 x 10" 3 (for 60% 
infrared flux accounted) and ^3 x 10~ 6 -5 x 10~ 3 (for 100% infrared flux 
accounted). The estimated ratio is smaller than the normal Galactic value (~ 
0.01; Spitzer, 1978), implying that a fraction of the shocked dust grains have 
been destroyed by sputtering. Given the lifetime of the dust grains destroyed 
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by sputtering for gas temperature ^ 10 6 K, r sput ~ 10 6 a(/zm)/n(cm~ 3 ) yr 
(where a is the grain radius; Dwek & Werner 1981), and adopting the gas 
density n ~2 cm~ 3 for the 3C 391 (CSSW04), we could essentially expect that 
some large grains (a <; 0.01/xm) which were initially present survive for the age 
of the remnant (~ 4x 10 3 yr) (CSSW04). 



4 Summary 

We have compared the morphology of SNR 3C 391 with the data of Chan- 
dra X-ray, SST mid-IR, and 1.5 GHz radio observations. The X-ray surface 
brightness is generally anti-correlative with the IR and radio brightness. The 
combination of the three band data clearly exhibits an multi-shell structure 
and a heart-shaped entire morphology for the remnant. Some detailed struc- 
tures are revealed: 

1. A thin brace-like shell detected at 24 /zm, which is projected outside the 
radio border confines the southern faint X-ray emission. The cusp is coincident 
with sn X-ray protrusion. 

2. The 24 /im knot on the SE boundary appears to be partly surrounded by 
soft X-rays. 

The mid-IR emission is probably dominated by dust grains, which are heated 
by the shocked hot gas and may have been partly destroyed by sputtering. 
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